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Abstract 

The relation for transition form factors of ij and rj' mesons is obtained by combining 
the exact nonperturbative QCD sum rule, following from the dispersive representation 
of axial anomaly, and quark-hadron duality. It is valid at all virtual photon momenta 
and allows one to express the transition form factors entirely in terms of meson decay 
constants. This relation is in a good agreement with experimental data. 
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1 Introduction 

Axial anomaly [Tl[2] is known to be a fundamental notion of nonperturbative QCD and 
hadronic physics. Usually it is considered for the case of real photons, however, the dispersive 
form of axial anomaly [3J can be considered for virtual photons also |1H6]. It leads to an 
exact anomaly sum rule (ASR) which does not have perturbative corrections due to Adler- 
Bardeen theorem [7] as well as nonperturbative QCD corrections due to 't Hooft consistency 
principle. Recently, this sum rule was applied to the analysis of pion transition form factors 
[8] which allowed to validate the interpolation formula for pion transition form factor [9l[T0] . 
This form factor attracted much attention because of unexpected and provocative data of 
BABAR collaboration [TT]. These data were followed by a stream of theoretical papers, 
that in particular questioned the QCD factorization [121113], provided the detailed analysis 
of perturbative and nonperturbative QCD corrections [HHTT] and suggested various model 
approaches jT8] - [22] . In our anomaly-based analysis [8] it was found that the BABAR puzzle 
may indicate an existence of small nonperturbative correction to continuum which must be 
compensated by significant correction to transition form factor. 

Recently, the BABAR collaboration extended the analysis and presented the data for t] 
and 7]' meson transition form factors |23j. These data motivated several recent papers [2H428] . 

In this work we analyze the r] and rj' transition form factors by means of generalized ASR 
to include the effects of meson mixing. Using the dispersive representation of axial anomaly 
the particular (octet) combination of meson transition form factors is expressed in terms of 
meson decay constants only. This expression is in a good agreement with experimental data 
in the whole range from real to highly virtual photons. 

2 Anomaly sum rule in octet channel 

In this section we briefly remind the dispersive representation for axial anomaly [3] (see 
also [2S] for a review) and derive anomaly sum rule for the octet channel of axial current. 
The VVA triangle graph correlator 



contains axial current J^s and two vector currents = {euU'-yf^u + Cdd'y^d + esS'j^s); k,q are 
momenta of photons. The octet component of axial current which is relevant for us explicitly 
reads: J^J = ^(M7a75« + d'^a'^bd - 257^755). 

It is convenient to write the tensor decomposition [301432] of correlator ([1]) in a form: 



where the coefficients Fj = Fj{k^, q^.p^] rn?), p = k + q, j = 1, . . . ,6 are the corresponding 
Lorentz invariant amplitudes constrained by current conservation and Bose symmetry. Note, 




(1) 




Fl Safiupk^ ~\~ F2 Safiupq^ 

+ F3 k^Eappak'^q" + Fi q^Eappak^q" 
+ F5 k^Eaupak'^q" + Fq qpEaupcTk''q" , 



(2) 
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that the latter includes the interchange fi ■H- u, k ■H- q in the tensor structures and fc^ f-> 
in the arguments of scalar functions Fj. 

In this paper we are interested in the case of one real (/c^ = 0) and one virtual photon 
(Q^ = —q^ > 0). Then for the invariant amplitude — Fq the anomaly sum rule (ASR) 
takes the form [S]: 

/•OO 1 

/ AUt;q',m')dt=—N,C'~''\ (3) 
where iVc = 3 is a number of colors and 

A3a = ^Im{Fs-Fe), (4) 

C<».^-L(4 + e3-2e;)^^. (6) 

The ASR ([3]) is an exact relation, i.e. it does not have perturbative corrections to the 
integral [7] as well as it does not have any nonperturbative corrections too (as it is expected 
from 't Hooft's principle). Another important property of this relation is that it holds for 
an arbitrary quark mass m and for any g^. 

Note, that one can also write out the similar relation with 

+ (7) 

for the singlet component of the axial current. However, in this case the absence of corrections 
is not guaranteed and one should explicitly take into account the gluonic anomaly. So we 
will not consider the singlet channel in this paper concentrating on the octet one. 

Saturating the l.h.s. of the 3-point correlation function ([1]) with the resonances and 
singling out their contributions to ASR ([3]) we get the sum of resonances with appropriate 
quantum numbers: 



f^F^ + f^iFrf/ + {other resonances) = 

poo T 

/ AUt:q\m^)dt = —N,C^'\ (8) 

Here the projections of the axial current j'f^ onto one-meson states M(= rj, rj') define the 
coupling (decay) constants flf. 

{0\jSmM{p))=tpJt,, (9) 
while the form factors Fm^ of the transition 77* — )■ M are defined by the matrix elements: 

d'xe''^M{p)\T{J,{x)MO)}\0) = e,,,^kPq''FM, • (10) 
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The relation ([HD is exact and expresses the global duality between hadrons and quarks. 
Nevertheless, to analyze the hadron properties one should additionally implement the local 
quark-hadron duality hypothesis. Taking into account large t] — rj' mixing, one can express 
the spectral function A^aiSjQ"^) in form of "two resonances+continuum" : 



A3a {s, Q') = vr/«5(s - mJ)F,^ {Q') + 
7if^,6{s - ml,)F,,, {Q') + - so). (11) 



Here Sq is a continuum threshold and A^^^ at one-loop level is 

.QCD _ 1 /.^x 

As it was shown in p3] there is no two-loop as corrections to this expression. 

Analyzing (11 II) . ( JT2l) one should note that the particles with nonzero two-photon decays 
cannot be included in the continuum as it vanishes at = 0, so they should be taken 
into account explicitly in the ASR. For heavy mesons the corresponding coupling constants 
should be suppressed [SIEH] at least as (m.^/m.res)^ which follows from the conservation of 
axial current jj^^ in the chiral limit (if only strong interaction is taken into account). That 
is why we restrict ourselves only to rj and rj' mesons. The ASR for the octet channel then 
reads: 

7r/«F,,(g2) + vr/« F,,,(g2) = ^ /° . (13) 

ZTTVO W + So 

Let us stress that this relation is correct for all due to the absence of the corrections 
to the Im{F'3 — Fq) [33] which allows to utilize the above expression for different . 

For real photons (Q^ = 0) the above expression coincides with the expression in [35], 
which was obtained from dispersive approach to axial anomaly in somewhat different way: 

7r/^«F,,(0) + vr/« F,,,(0) = (14) 

where 



FM,{^) = \r^V¥- (15) 

Here TM^2'y and tum are two-photon decay widths and masses of 77,77' mesons correspond- 
ingly. 

The equation (fT3|) allows us to fix the continuum threshold Sq by considering the limit 
— 00 where the QCD factorization [SSfl^ is applicable (contrary to ASR the exploration 

of this limit in generic QCD sum rules is obscured by possible corrections). The form factors 

at large [?,|38l|39] are: 

Q2pas ^ 2(c(«)/« + C(°)/°) (16) 
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Q'F^r^ = 2(CW/« + CW/J) /' (17) 

Jo ^ 

We take into account that in the hmit ^ oo the hght cone distribution amphtudes of 
both ri,T]' mesons are described by their asymptotical form [36|[37]: 0''*(x) = 6x(l — x). 
Then the ASR for the octet channel at large leads to: 

^0 = 4vr2((4)^ + (40' + 2x/2[/«/^° + /«/;]). (18) 

Substituting (fT8|) into ( IT3|) we express ASR in terms of meson decay constants fpj only, 
which is our main result: 



f^F,,{Q') + f^,F,,,{Q') = 

V 3 

4vr2 + QV((/,^)^ + if,')' + 2v^[/,«/° + /,V°.]) 

It is instructive to compare this formula with interpolation formulas for transition form 
factors of rj, rj' mesons proposed in [H]. The use of the proposed interpolation formulas 
leads to relation which is different from ( IT9l) but coincides with it in the limit of large and 
small Q^. Numerically, the difference is small at all (the maximal difference is 10% at 
~ 1 GeV) which proves that the interpolation formulas are good approximations. 

Let us now pass to apphcations of f|T9l) . 



(19) 



3 Mixing schemes and transition form factors 



Let us analyze the ASR (1191) for different mixing schemes. It is convenient to eliminate the 
dependence of the l.h.s. of f|T9|) on the scale of the coupling constants by dividing both sides 



of Eq. ( fT9l) by fs = \J [f^Y + WfiiY- Consider the form factors multiphed by 

^\f,F,,{Q') + f,,F,,,{Q')) = 



2 

/8 V 3 



4vr2 + QViif^y + Ul'Y + 2v^[/«/° + /J/J]) 
and introduce the matrix of decay constants in the usual way: 



— ( fri fri' 



We start from the one-angle mixing scheme, with 

fs cos fs sin t 
-fosinO focos9 



(20) 



F = ( 7o 7o ) ■ (21) 



p ^ ( /scos^ /gsin^ , 
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Figure 1: ASR for one- angle mixing scheme 
(pgj): e = -14° 
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Figure 2: ASR for one-angle mixing scheme 
([23]): e = -1Q° 



This mixing scheme was analyzed in many papers (see e.g. [12] and references therein), 
giving the values of mixing angle in the range 6 = —12° -. — 22°. 
In this case the ASR acquires a simple form: 

Q\F,,{Q') cose + F,,,{Q') sine) = ^^-^j^^^^, (23) 

where constant fg is defined by the anomaly sum rule at = f|T4|) : 

fs = ^-(,l^^cos9+,l^^sm9] . (24) 

Thus, (123|) and determine the mixing angle in terms of physical quantities (decay 
widths and transition form factors). 

The corresponding relation is plotted for different mixing angles in Fig. 1,2. The dots 
with error bars correspond to the l.h.s. of Eq. ( I23l) . where the form factors of 7],7]' mesons 
are taken from experimental data of CLEO |13] and BABAR [23] collaborations. The r.h.s. 
of Eq. fl20]) corresponds to the curve with the shaded stripe defined by the experimental 
uncertainties of meson decay widths. The slope of the straight line in the origins indicates 
the value of the anomaly sum rule at = 0. This sum rule was considered earlier [35j and 
is in a good agreement with the experimental values of two-photon decay widths of rj and t]' 
mesons. 

We observe that for 6 = —14° -. — 16° there is a reasonable agreement with the experi- 
mental data [23]H3]. The corresponding value of coupling constant is fs = 0.87^0.94/^, 
U = 0.13 GeV. 

The slight increase at large is related, in particular, by the tendency of rj' contribution 
to decrease at large Q^. Because of the negative mixing angle this leads to the behavior in the 
octet channel, qualitatively resembling the pion case, which produced the BABAR puzzle. 

• Let us now discuss the mixing schemes suggested and developed in |m - H6] . These 
schemes parametrize the decay constants f^j in terms of two mixing angles 6s, do'- 
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Figure 3: ASR for scheme 
1.17/^, /8 = 1.26/^,^0 = -9.2°, 
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Figure 4: ASR for scheme 
1.29A,/8 = 1.51^,0o = -2.4°,5, 



40 



fo = 
-23.8° 



F=( ^r'^^ {'''''^'V (25) 
\^ -/osm6'o /ocos6'o / 

Note, that this kind of matrices may appear when one considers the quark basis (see, 
e.g. [15]). For the parameters, suggested in gS] {fo = 1.17/^, /g = 1.26/^, = -9.2°, = 
-21.2°) the plot describing the ASR ([20]) is shown in Fig.3. 

For parameters, suggested in 06] {fo = 1.29/^, /§ = 1.51/^, = -2.4°, = -23.8°) the 
plot for ASR (1201) is shown in Fig. 4. One can see that for this particular mixing model the 
agreement of ASR with the experimental data is substantially worse. 

The developed approach can be used as an additional constraint on mixing parameters. 
It is useful to analyze the mixing schemes which take into account more than two mixing 
states; the detailed analysis of these schemes is quite lengthy and will be presented elsewhere. 



4 Discussion and Conclusions 

We generalize the rigorous nonperturbative QCD approach which relies on quark-hadron 
duality hypothesis but does not contain any free adjustable parameters to the case of rj and 
7]' mesons where mixing is crucially important. 

Combining the exact dispersive form of anomaly relation, quark-hadron duality hypoth- 
esis and asymptotic matching with QCD factorization we express the combination of i], rj 
meson transition form factors in terms of meson decay constants only (fTOl) . The obtained 
anomaly sum rule is valid in the whole kinematical region starting from = 0. This ASR 
is quite robust and can be used as a test for different sets of mixing parameters. Our anal- 
ysis shows that for a large number of mixing schemes the ASR is in a good agreement with 
experimental data (probably with the exception of parameters offered in [16]). 

At the same time, let us note that if one estimate the value of continuum threshold so 
from Eq. f|T8l) it appears to be quite small for all mixing schemes. This may reflect the 
contradiction of local quark-hadron duality from one side, and anomaly from the other side. 
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The possible resolution of such contradiction can be due to suppressed nonperturbative 
corrections to continuum similar to suggested earlier ^ for explanation of BABAR pion 
puzzle. It is significant that such corrections are essential for the value of sq while ASR 
ffTOj) is practically insensitive to them. This problem as well as the consideration the singlet 
channel of axial current and other mixing schemes requires a special investigation which is 
now in progress. 

We thank B. L. loffe, P. KroU, S. V. Mikhailov, T. N. Pham, A. V. Radyushkin for useful 
discussions and elucidating comments. This work was supported in part by RFBR (Grants 
09-02-00732, 09-02-01149, 11-02-01538,11-02-01454) and by fund from CRDF Project RUP2- 
2961-MO-09. 
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